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ABSTRACT 

We present the first detection of extragalactic submillimcter H2O maser in the 32f GHz 
transition toward the center of Circinus galaxy, the nearby Type 2 Seyfert using the Atacama 
Large Millimeter /Submillimeter Array. We find that Doppler features of the detected 321 GHz 
H2O maser straddle the systemic velocity of the galaxy as seen in the spectrum of the known 
22 GHz H2O maser in the galaxy. By comparing the velocities of the maser features in both 
transitions, it can be deduced that the 321 GHz maser occurs in a region similar to that of the 22 
GHz maser, where the sub-parsec-scale distribution of the 22 GHz maser was revealed by earlier 
very long baseline interferometry observations. The detected maser features remain unresolved 
at the synthesized beam of ~ 0".66 (~15 pc) and coincide with the 321 GHz continuum peak 
within small uncertainties. We also present a tentative detection of the highest velocity feature 
(redshifts up to ^635 km s -1 ) in the galaxy. If this high- velocity feature arises from a Keplerian 
rotating disk well established in this galaxy, it is located at a radius of ~ 0.018 pc (~1.2 x 10 5 
Schwarzschild radii), which might probe molecular material closest to the central engine. 

Subject headings: galaxies: active — galaxies: individual (Circinus) — galaxies: ISM — galaxies: nuclei 
— masers — submillimeter: galaxies 



1. INTRODUCTION 

The H2O masers in the 6i6~523 transition (rest 
frequency (f re st) = 22.23508 GHz) in active galac- 
tic nuclei (AGNs), what we call megamasers, have 
been detected in ~ 150 galaxies, most of which 
are detected exclusively toward narrow-line AGNs 
having a Seyfert 1.8-2. or LINER nucleus (e.g. 



spectra giving evidence for emission from edge-on 
rotating disks on scales of sub-parsecs as mea- 
sured in NGC 4258 with very long baseline in- 
terferometry (VLBI) observati ons at milliarcsec- 



ond angular resolution (e.g., Mivoshi et al. 19951 : 
Hagiwara etail l200ll iKuo et al.l l201lh . Accord 
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Kondratko et al.l 120061 lHagiwara et al 
Hagiwarall2007l ). Some of these megamasers show 
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ing to theoretical studies, the 22 GHz H2O maser 
is inverted in wider range of physical conditions 
at kinetic temperature of T k = 200-2000 K, 
hydrogen densities of n(H2) = 10 8 -10 10 cm~ 3 , 
and H 2 densities of n(H 2 0) = 10 3 -10 5 cm" 3 , 
which largely excites submillimeter H2O maser 
transitions such as the 321 GHz (1029-9.3 6) and 



325 GHz (5i5-4 22 ) transitions (IDeguchil 1977 



i Neufeld fc Melnick|[l990tlNeufeld fc Melnicklll991 
lYates et alJll997|) . These studies have shown that 
the 321 GHz maser is even strongly excited un- 
der more restricted physical conditions of Tk > 
1000 K. The upper state energy of the 321 GHz 
transition is higher than the 22 GHz transition 
(£? u /fc=1862 K for the 321 GHz transition and 
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E u /k=6U K for the 22 GHz transition). 

The first discovery of extragal actic sub-millimete r 



H2O masers was achieved by iHumphrevs et al 



(|2005l ). in which the detection of H2O masers in 
the 183 GHz (3i 3 -2 20 ) and 439 GHz (6 43 -5 5 o) 
transitions was reported in the nuclear region of 
NGC3079, the Type 2 Seyfert/LINER galaxy, al- 
though the detection of the 439 GHz maser was 
tentative. The H2O emission in the 183 GHz 
transition was also detec ted toward the luminou s 
infrared galaxy Arp220 (jCernicharo et al" 2006). 
However, the detected emission could be thermal 
emission. 

Circinus galaxy at a distance of 4.2 ± 0.8 
Mpc is a nearby spiral l ying b ehind the Galac- 
tic plane! Freeman et"al1 Il977l) . exhibiting evi- 
dence of the co-existence of starbursts occur- 
ring within a nuclear mo lecular ring in CO emis- 
sion (jCurran et al.l 120081 ). The galaxy hosts a 
Type 2 Seyfert nucleus having a luminous (~ 
10 42 erg s -1 ) hard X- ray source obscu red by 
Compton-thick medium Matt et al. ( 19961 ). The 
adopted systemic velocity of the galaxy is 433 



km s 1 (de Vaucouleurs et al.l 1 1991; Braatz et al 



2003). (The optical velocity definition is adopted 
throughout this article and the velocities are cal- 
culated with respect to the local standard of rest 
(LSR).) 

Bright 22 GHz H2O maser emission in the 
galaxy was first reported bv lGardner fc Whiteoak 
(| 19820 and the known highly variable flux den- 
sity o f the maser is in the r a nge of ~5 to ~40J 



200 



(e.g., Greenhill et al.l 1997 ; iBraatz et al 
McCallum et al.l 120051 ). The velocity range of 
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highly Dopplcr-shifted (high-velocity) maser fea- 
tures detected to date is ~50 to ^900 km s _1 
(jGreenhill et al.l l2003bl) . Follow-up VLBI obser- 
vations of the 22 GHz maser at the milliarcsecond 
resolution revealed that the maser traces a warped 
edge-on disk that implies a sub-Keplerian rotation 
with a radius of ~0. 1-0.4 pc at the maximum rota- 
tion velocity of ^270 km s" 1 around a nucleus and 
outflowing components distributed near the disk 



with t he velocity of ^160 km s 1 (jGreenhill et al 



l2003ah . The Circinus galaxy, already known to 
show the strongest flux density of the 22 GHz 
H2O maser among the H2O megamasers, could 
show the 321 GHz H2O emission as H2O masers in 
both the 22 and 321 GHz transitions have been de- 



tected in the same Galactic sources (Menten et al 



1990t iMenten fc Melnicflll99lh . 

In this work, we present the observation of 

the 321 GHz H2O maser emission in Circinus 

galaxy, obtained using Atacama Large Millimc- 

ter/submillimctcr Array (ALMA) in the Cycle 

program. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

We observed H2O emission in the IO29-936 
transition (^ res t=321.226 GHz) toward the Circi- 
nus galaxy on 2012 June 3 using the ALMA 
as part of a search for extragalactic submil- 
limeter H2O maser. Observations lasted for 49 
minutes (approximately 25 minutes on-source) 
with 18 antennas. The tracking center of the 
galaxy was a(J2000) = 14 h 13 m 09 s .906, (5(J2000) 
= -65°20'20".468. These observations were made 
with a single dual polarization spectral window 
(1.875 GHz bandwidth), divided into 3840 spec- 
tral channels, yielding spectral resolutions of 488.3 
kHz or 0.458 km s _1 at the observed frequency of 
319.4 GHz. A resultant total velocity coverage is 
- 1760 km s _1 . The 1.875 GHz spectral window 
was centered near the galaxy's systemic velocity of 
433 km s -1 . Phase- referencing observations were 
conducted, by switching to a phase-referencing 
source, J 1329-5608, the position of which is « 
9 away from the target source. Data calibra- 
tion was performed using the Common Astron- 
omy Software Applications (CASA). Amplitude 
calibration was performed using observations of 
Titan, and the bandpass correction was made with 
3C 279. Flux-density errors of 10 % are conserva- 
tively adopted, based on the capabilities of Cycle 
observations in Band 7. Imaging was performed 
using CASA in natural weighting, for which the 
synthesized beam size was 0".66 x 0".51 (P.A.= 
-17°). The image analysis was performed using 
both the CASA and the Astronomical Image Pro- 
cessing Software. After the phase and amplitude 
calibrations, the 321 GHz continuum emission of 
the Circinus was subtracted from the spectral line 
visibilities using line-free channels prior to the 
imaging and CLEAN deconvolution of line emis- 
sion. The line emission in the Circinus galaxy was 
thus separated out from the continuum emission. 
The rms noise level of spectral line maps in natu- 
ral weight was ~9— 11 mJy beam -1 per spectral 
channel, depending on each frequency channel. 
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The rms of the continuum map was ~0.8 mJy 
beam -1 . 

3. RESULTS 

The 321 GHz H2O maser in the Circinus galaxy 
was searched within the nominal field of view of 
18", centered on the phase-tracking center of the 
galaxy and in the velocity range of V^lsr=— 430 to 
+1290 km s , covering the total velocity range 
of known 22 GHz H2O m aser measured in ear- 
lier si n gle-dish obse r vation s ( Gardner & Whitcoak 



the systemic velocity, consistent with the fact 
that earlier observations detected no distinct 
systemic fea ture of the 22 GHz maser in the 



19821: iNakai et all Il995t IGreenhill et al 



Braatz et all 120031 : iMcCallum et al.l 120051 ). We 
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detected the 321 GHz H2O line emission, the 
peak flux density of which is 131 mJy beam -1 at 
Vlsr=531.1 km s -1 with a signal-to-noise ratio 
(S/N) of >10 (Figure [J) and the total flux density 
is ^350 mJy. The total integrated intensity esti- 
mated from all the detected H2O emission is ~14 
Jy km s . The apparent luminosity of the emis- 
sion is ~ 5 Lq, assuming isotropic radiation of the 
emission. The detected H2O emission can be di- 
vided into two groups, one of which is blueshifted 
emission of Flsr='~255-'^275 km s -1 peaking 
at ^lsr=274.4 km s -1 and the other is red- 
shifted emission lying from Vlsr=~470-~670 
km s -1 . The line shapes of these Doppler-shiftcd 
features are narrow, similar to the H2O maser 
features seen in the 22 GHz single-dish spectra 
(Figure [2]) . The redshifted emission is composed 
of broad features, narrow-line features peaking 
at Vlsb=477.1, 511.4, 530.6, 531.1, 541.6, 547.1, 
559.0, 587.8, 598.8, 602.4, 618.0 km s -1 , and a 
minor peak at V"lsr= 668.7 km s -1 (denoted by 
an arrow in Figure [2]). Figure displays a highly 
redshifted feature tentatively detected in ~ 3<r 
level centered at Vlsr=1069.1 km s -1 . One may 
suggest a possible detection of a feature peaking 
at V^lsr= 1129.7 km s -1 (< 3a), which must be 
confirmed in further measurements. 

In Figure [5J we made a comparison of the 
321 GHz H2O spectra with 22 GHz spectra 
obtained on 2012 September 7 using the Deep 
Space Network 70 m antenna at Tidbinbilla. The 
known strong variability of the 22 GHz maser 
in this galaxy makes it difficult to compare pre- 
cisely the velocities of each feature between the 
22 GHz and 321 GHz lines. In our 321 GHz 
observation, no feature is detected at or near 



galax y (e.g., IGreenhill et al.l l2003bt iBraatz et al 



l2003h . It should be also remarked that the asym- 
metry between the blueshifted and redshifted 
emission is observed as in the cases of other 
H2O megamasers such as NGC 4258: Redshifted 
features are more numerous and are significantly 
more intense than the blueshifted features (e.g., 
Maoz fe McKee1ll998h . Thus, the 321 GHz emis- 



sion in the galaxy largely overlap Doppler fea- 
tures appearing in the 22 GHz maser spectra ob- 
tained in the earlier single-dish observations. The 
Gaussian-fitted position (relative to the phase- 
referencing source) of the redshifted peak emission 
at 531.6 km s -1 is a(2000):14 h 13 m 09?95±0.01, 
5(2000): -65°20'20".92±0.02. The position of the 
blueshifted peak at 274.2 km s -1 is 
a(2000):14 h 13 m 09!96±0.03, <5(2000): -65°20'20".89 
± 0.03. The relative positions of other redshifted 
and blueshifted features coincide with these po- 
sitions within uncertainties of fitting errors and 
these features remain unresolved at the angular 
resolution of ~ 0".66, or 14.6 pc. 
The 321 GHz sub-millimeter continuum was also 
detected in the center of the galaxy with an 
S/N of > 50 (Figure |4}, the peak position of 
which is a(2000):14 h 13 m 09.96±0.01, 5(2000): 
65°20'20.87"±0.01. The detected nuclear contin- 
uum with the total flux density of 55 mJy was 
partly resolved in the northeast direction accom- 
panying an elongated substructure in the north 
(Figured]). According to the ALMA Cycle Capa- 
bilities, the positional errors are ~ 10% of the syn- 
thesize beam, 0". 66/10 or 0".066 in our observa- 
tion, which seems to be the most dominant source 
of an error. Approximate relative positional er- 
rors between the H2O emission and the contin- 
uum peak, represented by the sy nthesized beam 
size d ivided by 2 x S/N (e.g., iHagiwara et al. 
2003bl ) are ~ 0".07 or 1.5 pc. Consequently, the 



relative positions between the maser and contin- 
uum emission peak are consistent within these 
errors. Note that the positions of the blueshifted 
and redshifted maser peaks at Flsr=274.2 and 
531.6 km s -1 are offset from the optical nu- 
cleus at a(2000):14 h 13 m 09.950±0.335, (5(2000): - 
65°20'21.20"±0.25 (primary from NED) by -0".28 
to -0".31 in declination, although these offsets arc 
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almost within errors. 

The resolved 321 GHz continuum emission shows 
an elongated structure (Figure |4j, the direction 
of which is similar to that of a molecular or 
atomic gas disk at an inclination of ~ 60°-70° 
(ICurran et alj|2008ft . 



4. DISCUSSION 

4.1. Origin of the 321 GHz H 2 and Sub- 
millimeter Continuum Emission 

The synthesized beam size of our observation 
was too large to determine if the detected 321 GHz 
H 2 emission is thermal emission or maser, from 
the value of brightness temperature. Given that 
the 321 GHz line features are confined within a 
diameter of <14 pc and the line shapes are nar- 
row (<0.45 km s -1 ), it is most likely that the 
H 2 emission detected in the Circinus galaxy is 
due to maser amplification. 

To investigate the nature of the 321 GHz contin- 
uum, more information such as a spectral index 
would be required. We can speculate about the 
possibility of thermal emission from dust as the 
continuum shows an extended part in the general 
direction of a galactic-scale disk in a dusty envi- 
ronment, while a central nuclear continuum could 
contain the contribution of free-free emission. 

Since the peak position of the continuum co- 
incides with that of the maser, the maser could 
be enhanced through the thick dust lane in our 
line of sight. The optical studies revealed that the 
central 15 pc of the galaxy is dominated by the 
AGN luminosi ty, with star formatio n contribut- 
ing only 2% (|Sosa-Brito et al.l l200lh . The 321 
GHz maser in the Circinus galaxy could be as- 
sociated with AGN activity. Because of the rel- 
atively smaller maser luminosity, we still cannot 
still rule out that the maser originates in star- 
forming activity. The distinct velocity gradient of 
the 321 GHz maser is not detected throughout our 
analysis of the first moment map produced from 
the spectral-line cube containing features lying at 
Vlsr=256-670 km s^ 1 , suggesting that the maser 
emission is neither spatially nor kinematically re- 
solved in our observation on scales of 10 pc, while, 
according to the earlier VLBI observations, the 22 
GHz mas er in the galaxy was resolved on scales 
of pal pc (jGreenhill et alJl2003bl : iMcCallum et al. 
2009), implying that a better angular resolution 



is required, at least, better by one order of mag- 
nitude than that of our observation. The future 
VLBI mapping of the 321 GHz maser at milliarc- 
second resolution should be able to clarify sub- 
parsec-scale kinematics of the galaxy by measur- 
ing distributions of each maser spot. 

4.2. Comparison of 321 GHz H 2 Emission 
with 22 GHz H 2 Maser 

In Figure 2, one can compare the spectra of the 
321 GHz maser emission with the 22 GHz maser 
emission, both of which were observed at approxi- 
mately the same time in 2012 June and September. 
The velocity ranges of masers at both transitions 
overlap to some extent and some of the prominent 
maser features peaked at or near the same veloci- 
ties (see Figure 2). However, the velocity span of 
the 321 GHz maser is smaller than that of the 22 
GHz maser due to the spectrum haying a broader 



velocity coverage ([Greenhill et al.l l2003bl ); high- 



velocity features of the 22 GHz transition span 
from ~50 km s _1 to ^900 km s _1 having the 
most Doppler-shifted velocity of 4 60 km s" 1 wi th 



respec t to the systemic velocity (jGreenhill et al 



2003bf ). while those of the 321 GHz transition 
from ^270 km s _1 to ^670 km s _1 in our obser- 
vation, excepting the minor detection. This might 
be explained by the fact that the 321 GHz maser 
is present in more physically restricted regions 
in higher temperature and density as predicted 
in theoretical stu dies ( Neufeld &; Melnick 1990; 
Yates et al.lll997l ). which should be examined by 



resolving the spatial distribution of the maser at 
higher angular resolution. However, with the com- 
parable sensitivity of the earlier 22 GHz maser 
spectra, the total velocity range of the 321 GHz 
maser emission might be as broad as that of the 
22 GHz maser. In the case of the Cepheus A 
star-forming regions, the 321 GHz and 22 GHz 
masers were observed within the region of ~ 1", 
in which spatial coincidence between the 321 GHz 
and 22 GHz masers was not found and the 321 
GHz maser is believed to be tracing warm er gas 



than the 22 GHz maser (jPatel et al.l 120071 ). We 



note that the 321 GHz maser in evolved stars is 
highly variab le, compared with 22 GHz and 325 
GHz masers ( Yates et al 1 ll996h . Thus, the study 
of the intensity variability of the 321 GHz maser 
in the galaxy is of great interest. 
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From Figure 2, it is difficult to find one-to-one 
correspondence of each velocity feature between 
the two transitions. However, the overall velocity 
spread of the 321 GHz maser features looks sim- 
ilar to that of 22 GHz masers that are believed 
to trace a masing disk rotating around a central 
engine having 1.7 x 10 6 M Q enclosed mass and 
non-disk maser components from bipolar outflows 
lying within about 1 pc from the edge-on disk 
(jGreenhill et al .1120 03a') . At the angular resolution 
on our map, it is not possible to distinguish only 
from their velocities whether or not each of the 
detected 321 GHz maser features is coming from 
the disk or non-disk outflows. 

4.3. Possible Detection of High-velocity 
Gas in the Circinus 

The known highest velocity maser feature in 
Circinus galaxy is redshif ted by ^460 k m s _1 from 



the systemic velocity (jGreenhill et al.l l2003bl ). 



while the blueshifted counterpart was not detected 
in our 321 GHz spectra. Assuming that the 1069 
km s _1 feature (redshifts up to ~ 635 km s _1 ) 
is real, it is located at a radius of ~ 0.018 pc 
(~ 1.2 x 10 5 Schwarzschild radii for a 1.7 x 10 6 
Mq black hole) from the central engine, by calcu- 
lation based on a Keplerian (v oc r~° ) rotation 
disk model in Greenhill et al. (2003a, 2003b). 
This high-velocity feature could probe the inner 
most part of a molecular gas disk in the galaxy. 
The submillimeter maser is potentially a powerful 
tool for exploring regions which cannot be probed 
with the 22 GHz maser. 



The interpretation of the origin of the H2O maser 
in the Circinus galaxy is not straightforward. Nev- 
ertheless, our discovery of the 321 GHz H2O maser 
in Circinus galaxy will be useful as a basis for 
future detailed studies of circumnuclear gas in 
AGNs. The detected 321 GHz and well-studied 22 
GHz traditions are from ortho-H20. Using these 
two different maser transitions will enable us to 
place new constrai nts on radiative transfer models 
( Humphreys! 2007 ) to determine physical condi- 
tions of molecular gas in the inner sub-parsecs of 
AGNs. Moreover, the H 2 lines at 22 and 321 
GHz show higher atmospheric transmission than 
those at 183 and 325 GHz, which is advantage for 
interferometric observations. Ultimately, we hope 
to measure the angular distribution of submillime- 



ter H2O masers including those of the 321 GHz or 
other transitions on the basis of currently on-going 
efforts in developments of VLBI at submillimeter 
wavelengths. 

This Letter makes use of the following ALMA 
data: 2011.0.00121. S. ALMA is a partnership of 
ESO (representing its member states), NSF (USA) 
and NINS (Japan), together with NRC (Canada) 
and NSC and ASIAA (Taiwan), in cooperation 
with the Republic of Chile. The Joint ALMA Ob- 
servatory is operated by ESO, AUI/NRAO and 
NAOJ. The Tidbinbilla 70m telescope is part 
of the NASA Deep Space Network and is oper- 
ated by CSIRO. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under 
contract with the National Aeronautics and Space 
Administration. Finally, we thank the anonymous 
referee for helpful suggestions that significantly 
improved the article. 
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Fig. 1. — Spectrum of the 321 GHz H2O maser from the center of the Circinus galaxy at the spectral 
resolution of 488.3 kHz or 0.458 km s _1 , obtained using 18 antennas on 2012 June 3 in the Cycle ALMA 
program. The total velocity range covered in the spectrum is Flsr= _ 300. 0-1250.0 km s _1 , outside of 
which no significant emission was detected. The LSR systemic velocity, 433 km s _1 is denoted by an arrow. 
Amplitude scale is in mJy beam . 
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Fig. 2. — Comparison of the H2O spectra between lines in the 321 GHz and 22 GHz transition toward the 
center of the Circinus galaxy. The 22 GHz spectra were obtained with the Deep Space Network (DSN) 70 
m antenna located at Tidbinbilla (S. Horiuchi et al., in preparation). The observations of the Tidbinbilla 70 
m were made on 2012 September 7 at the spectral resolution of 31.25 kHz or 0.42 km s _1 with an rms of 
^2.5 mJy and the 321 GHz spectra were on 2012 June 3. Minor features (Vlsr=206.4, 370.9, and 668.7 km 
s _1 ) are denoted by arrows. The flux density scales of the 22 GHz spectra shown in the figures are rather 
tentative and require detailed analysis (S. Horiuchi et al., in preparation). Upper left: blueshifted velocity 
range (200-400 km s _1 ) of the 321 GHz maser spectrum obtained from Figur^TJ Lower left: spectrum of 
the blueshifted maser obtained at Tidbinbilla. Upper right: redshifted velocity range (450-725 km s _1 ) of 
the spectrum from FigureflJ Lower right: spectrum of the redshifted maser obtained at Tidbinbilla. 
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Fig. 3. — Redshifted 321 GHz maser spectrum: a highly redshifted feature tentatively detected (~ 3<r ) 
peaking at Flsr= 1069 km s _1 and a possible detection of a feature at V^sr— 1129.7 km s _1 . Upper: 
amplitude of the spectrum between Vlsr= 1050 and 1150 km s _1 . Lower: corresponding phase. The phases 
of the 1069 km s _1 feature are indicated by a circle. 
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Fig. 4. — Integrated intensity image (contours) of the 321 GHz H2O maser overlaid the 321 GHz continuum 
image (gray scale) observed with the ALMA. The peak integrated intensity of the H2O maser and the peak 
flux density of the continuum are 13.6 Jy km s -1 beam -1 and 40.6 mJy beam -1 , respectively. Contour 
levels and gray scales (scales are from 2.7 mJy (~3<r) to 30 mJy) are denoted in the image. The synthesized 
beam (FWHM) is plotted in the left bottom corner. 
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